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ABSTRACT: For r-process nucleosynthesis the beta decay rates for a number of neutron- 
rich intermediate heavy nuclei are calculated. The model for the beta strength function is 



' able to reproduce the observed half lives quite well. 
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More than half of the nuclei heavier than iron are produced through r-process nucleosyn- 
O I thesis [1]. Though one knows that high neutron density and fast neutron capture time 

scales compared to beta decay are required for this process, astrophysical sites for them 
■ are not clearly identified. As the r-process path often reaches highly neutron-rich nuclei, 

collecting experimental information about them sometimes becomes prohibitively difficult. 
On the other hand near the neutron shell magic numbers depending on the shell gap [2] and 



^ • beta decay half lives the chain often has to wait long enough before proceeding to the next 



heavier nucleus [3]. These are called waiting point nuclei and and can be identified by their 
observed relative abundances. With the lack of experimental information, one often has to 
fall back on the predictions of theory for these neutron-rich nuclei. In this letter we consider 
such nuclei near = 82 magic shell for 115 < A < 140. The shell closure at A = 82 for 
n-rich nuclei in the r-process path is receiving a lot of attention lately and realistic models 
for beta decay half-lives and rates in this region are badly needed. We construct here a 
model for the beta decay strength function based on the experience with lighter nuclei in 
the fp shell. We deal with nuclei for which a few lowlying log ft's are known, to test the 
model and then also apply it to cases where not much is known experimentally beyond the 
Q-value and half life. The actual beta decay rates are also calculated for stellar density and 
temperatures which are thought to be typical for the r-process. 

Calculation of beta decay and electron capture rates have been pursued vigorously during 
the last twenty years for pre-supernova evolution and gravitational collapse stage for core 
collapse supernovae [4]. Whereas for the sd-shell nuclei the rates are given by Oda et al. 
[5], for the fp-shell nuclei for quite sometime the works of Fuller, Fowler and Newman 
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[6] were used as the standard refrence. Aufdcrheide et al. [7] indicated the usefulness 
of interacting shell model for the rates. They also calculated approximate abundances in 
stellar conditions appropriate for the pre-supernova and collapse stage and identified the 
nuclei which arc the most important in the network for each density-temperature grid point. 
On the other hand statistical methods for the beta decay strength distribution were used 
to calculate the rates for fp-shell nuclei [8]. This was followed by large space and detailed 
shell model calculations of the allowed beta decay strength to calculate the weak interaction 
rates for nuclei with A < 65 [4] [10]. The effect of the new shell model rates on the late 
stage of evolution of massive stars has also been seen [11]. Recently Pruct and Fuller 
[16] used FFN-type ideas supplemented by experimental information wherever available 
to calculate the rates for nuclei relevant for supernova collapse. However for much lower 
densities and temperatures appropriate for s- and r-process nucleosynthesis the work of 
Takahashi and Yokoi [17] is still being widely used. As large number of nuclei arc involved 
in the r-process path the estimates of half lives use a combination of global mass models 
and the quasi-particle random-phase approximation (QRPA). These are FRDM/QRPA [18] 
and ETSI/QRPA [19]. Self-consistent Hartree-Fock-Bogoliubov plus QRPA model have also 
been used [20]. 

The beta decay rate for a nucleus in stellar conditions at temperaure T is given by 
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where Ei is the energy of the state of the mother nucleus and Z is the partition function of 
the nucleus, 'j' sums over states of the daughter nucleus to which transitions are allowed. 
The rate from the parent nuclear state 'i' to the daughter nuclear state 'j' is given by 



fij is the phase-space factor for beta decay with an electron gas outside [8]. For the 
Coulomb correction factor we use the Schenter-Vogel expressions [9] [8]. The allowed ft 
values have two components 
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For calculation of half lives the decay is from the ground state, i.e., there is only one 
mother state with the Boltzmann factor equal to 1 in eq (1). In eq (3) the first part is 

the contribution from the allowed Gamow- Teller operator and the second from the allowed 
Fermi operator. Normally the GT strength distributes among three different types of final 
states of the daughter : 

(a) a discrete set of low-lying states for which the log ft's arc known experimentally 

(b) another set of discrete states above them for which the strengths are not known 
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(c) a part of the GT giant resonance. 

Usually one sees only the tail of the giant resonance as the Isobaric Analogue State (IAS) 
for P~ transitions is pushed up by the Coulomb interaction. In our treatment we include in 
A of eq (1) as well as in the half life the states of type (a), in particular all the lowlying states 
connected by allowed transitions wherever known experimentally. For the contribution from 
(c), i.e., the giant resonance region, we adopt the procedure detailed below. 

The Fermi resonance is easy to take care of. In the absence of Coulomb interaction, the 
total Fermi strength goes to IAS. Coulomb interaction spreads this strength into a sharp 
resonance around the IAS and we take the resonance width as ar = 0.157Zyl-V3 [8] [12]. As 
the spreading width is small this resonance cannot be reached by Q-value and so Fermi makes 
very little contribution to half lives and rates. But the situation for GT is quite different. As 
the spin-isospin dependent nuclear Hamiltonian does not commute with the GT operator, 
the GT resonance is a broad one. Detailed shell model results [13] for fp shell nuclei and 
use of spectral distribution theory [14] indicate that the GT strength distribution to a good 
approximation can be taken as a Gaussian for the initial state somewhat high in excitation 
(a few MeV) to be in the chaotic regime. This point has been studied at length in the 
last few years. A good understanding of the Gaussian forms follows by using Hamiltonians 
from the two body random matrix ensembles which model chaos in many particle quantum 
systems very well [21]. Shell model results for lighter nuclei show that even with specific 
two body interactions the Gaussian results persist. However for the ground state region 
there is some departure. In this letter we do not take into account such departure and take 
the GT resonance for each nucleus given by a Gaussian. Hence the problem comes down to 
fixing the centroid and width of the resonance. For the GT centroid Bertsch and Esbensen 
[15] using the Tamm-Dancoff Approximation give the following expression [16] 

Egt- = ElAS + ^Es.o. + 2[k^rSGT- - {N - Z)kr] (4) 

where ^Es.o. is the contribution coming from the spin-orbit force. For AE'^.o. we use 
the orbit averaged value of 3.0 MeV [15]. The last two terms come from the spin-isospin 
and isospin dependent nuclear forces. Here for the sum rule strength Sqj-- we use the 
expression 3(A'^ — Z). Following Pruet and Fuller [16] we use the values kr = 28.5/^4 and 
kuT = 23/A. For the overall normalisation of the sum rule Gaussian we use the quenched 
3(iV — Z) value minus the experimental strength already included in (a). For the quenching 
factor for the beta decay we use the value 0.6 [8]. Finally to include the contribution from 
(b), i.e., the discrete states below the GT resonance region, we artificially broaden the 
resonance by using the width larger than the normal giant resonance widths. This is open 
to improvement. Perhaps partitioning of the space into subspaces and taking into account 
a number of Gaussian strengths instead of a single Gaussian is a better way of handling 
this. But then the width of each Gaussian either has to be calculated theoretically or used 
as a parameter. We first use this model to calculate half lives of nuclei around the N = 82 
shell closure. We include in this letter nuclei in the range 115 < ^ < 140 and also use the 
selection criterion that the Q-value greater than 5 MeV for the statistical treatment of the 
high lying states. We separate the nuclei into two groups - one set with known experimental 



3 



ft values and the other with no knowledge about low lying ft's. 

We also mention here that the spectral distribution result for the strength distribution gives 
it as a bivariate Gaussian in both the initial and final energies. However for the r-process 
nuclei the stellar temperature (in MeV) is quite low compared to the energies of the excited 
states of the mother. As a result the decay takes place only from the ground state of the 
mother and hence the full bivariate GT distribution does not come into play. In Table 1 
we present the calculated half lives for 13 nuclei compared to their experimental values. 
These are obtained by taking the width of the GT giant resonance as a free parameter 
and minimising the summed [log^Ty^/ry^^)]'^ . The calculated values show reasonably good 

agreement with the observed values and S[log(r{^^2/'^i/2*)]^/-^ ^-^^ where 'N' stands for 
the number of nuclei. The optimised width is 5.0 MeV. 

Table 2 shows a similar comparison for nuclei for which no lowlying log ft's arc known. 
Here the minimised T.[log{T'^'^yT''^jf)]'^/N is 0.07 for the width 7.7 McV. The fact that the 
width for this case is larger than the earlier one is expected as here the resonance has to be 
extended to the ground state region. 

In Table 3 we calculate the beta decay rates for ^^^Sn and ^^^In for the purpose of illus- 
tration over a range of density and temperature. The chemical potential of the electrons 
changes from 1.32 MeV for the density 10^ g/cc to zero for 10"^ g/cc for the temperature 
5 X 10^ K. At higher temperatures the chemical potential is less. For 3 x 10^ K it is only 0.062 
MeV for density 10^ g/cc. We have considered the density range far beyond the r-process to 
examine the dependence on density. In this range this dependence is found to be very mild 
but the temperature dependence is seen to be a bit stronger for ^^^Sn. Of course the rates 
can be calculated for other nuclei and extended to other densities and temperatures. As the 
values of kT are small compared to the excitation energies of even the first excited state of 
the mother nucleus we have given only the ground state contribution. For similar reasons 
the rates coming from 'back resonances' of some special excited states with large overlap 
with the daughter ground state region are also unimportant here and hence are neglected. 
In future we plan to extend these calculations to more neutron-rich nuclei for which very 
little experimental information is available. We shall also consider heavier nuclei. We thank 
A. Ray and M.M. Sharma for encouragement. 
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Mother Daughter Q- value No of lowlying ri/2 ^1/2 
nucleus nucleus log ft's taken Exp Calc 

(MeV) (sec) (sec) 



1. Ifigg i38xe 7.820 3 6.49 12.99 

2. |fl84 i37xe 5.880 2 24.5 33.1 

3. g%3 ^^^Xe 6.930 3 83.4 69.5 

4. ifSbga ^^^Te 8.420 2 10.43 7.72 

5. sfSngg ^3^Sb 7.830 4 1.44 1.19 

6. IfSbgi ^^^Te 5.290 2 167.4 141.3 

7. livings ^^^Sn 6.746 2 0.282 0.280 

8. Iflngi i30gj^ 10.250 3 0.32 0.25 

9. Iflujg i28g^ 898Q 3 2.30 

10. IfCdyy 125j^ 7;lqq 3 Q 55 2.12 

11. i?°Ag73 i20c^ 8200 4 1.23 0.987 

12. i^^Agyi ii^Cd 7.060 3 3.76 7.23 

13. IfRh^i ^^^Pd 8.900 3 0.68 0.49 



Table 1. Half- lives of nuclei for which log ft values are available and were used. 
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Mother Daughter Q- value ri/2 t-i/2 
nucleus nucleus Exp Calc 

(MeV) (sec) (sec) 



1. lfTe^5 ^^^I 6.370 1.4 1.86 



2. |fTe84 i37j 6940 2.49 1.41 



3. M^Sbgs ^^*^Te 9.300 0.201 0.208 



4. iflngg i^^s^ 0.82 0.507 



5. ifRh72 ^^■^Pd 7.000 0.44 1.35 



Table 2. Half-lives of nuclei for which no log ft values are available. 



7 



NUCLEUS DENSITY 

( g/cw? ) 



TEMPERATURE 
Tg (in 10^ K) 



0.5 1.0 2.0 3.0 



133 



Sn 



10« 
10^ 
10^ 
10^ 
10^ 
10=^ 



0.543 
0.581 
0.584 
0.584 
0.584 
0.584 



0.542 
0.558 
0.559 
0.560 
0.560 
0.560 



0.493 
0.498 
0.498 
0.498 
0.498 
0.498 



0.450 
0.452 
0.452 
0.452 
0.452 
0.452 



108 


1.70 


1.75 


1.59 


1.46 


10^ 


1.93 


1.83 


1.61 


1.47 


10^ 


1.95 


1.83 


1.61 


1.47 


10^ 


1.95 


1.84 


1.61 


1.47 


10^ 


1.95 


1.84 


1.61 


1.47 


10^ 


1.95 


1.84 


1.61 


1.47 



Table 3. Decay rates in s ^ for ^'^'^Sn and ^'^^In. 
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